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’ ThF MAGNETIC FIELD OF MERCURY 

A.iSTriAC 7 

Direct observations of the magnetic field of Mercury and its 
inngnotosphere, formed by the interaction of the solar wind, weiv 
performed twice by the Mariner 10 spacecraft, in March of 1974 and 
aguin in March 1975, From these data, it is clear that there cxis..*> 
an intrinsic uagneti''. field of the planet, sufficiently strong at 
present to deflect the solar wind flow around the planet and to fora 
a detached bow shock wave in the super Alfvenic solar wind. 

Four methods have been used to analyse the magnetic field data and 
uerive quantitative values for the description of the planetary field: 

1. Comparison of bow shock and magnetopause relative positions 
at Mercury to those at Earth, 

2. Direct spherical harmonic analysis of the data, 

3. Modeling of the magnetosphere by an image dipole and infinite 2-D 
current sheet in addition to the planetary field, 

4. Scaling of a mathematical model for the terrestrial 
magnetosphere. 

/ 02 1 2 12 

lire results obtained yield dipole moments^/ (g^ ) +(H^ ) "Kfh ) 

22 3 

ranging from 2.4 to 5,1x10 gauss-cm , with the lower values associ-ced 
with certain models using partial quadrupole (g^ 0 ) and octupole (g^ ) 
terns to improve the least squares fitting of models to observations, 
uocause the data set is incomplete, in the mathematical sense, no 
unique representation of the planetary field multipolar represent.** *on 
can be derived by method (2). The use of only 1 of the 5 quadrupole 
moment terms and 1 of the 8 octupole moment terms corresponds to a 
displacement of the dipole along its axis. These terms used in 


mvw'.oub (3) and (4), yield equivalent offsets of the dipole by 

approximately 0.2 Rj,.. The selection of only those higher ordc : torus 

possessing axial symmetry cannot be justified. Thus, the 

arge offset may reflect the limitations of the models used to 

represent the external current systems. Bccpu re of the relatively 

snort radial excursion of the data, the g 3 and terms can axso 

be spatially aliased with the g^° term. 

Analyses by method (2) of subsets of data from the third encounter, 

t«ken near closest approach, yield a convergent series of dipole iaOtitCUC 

values which are believed to best represent the intrinsic planetary 

3 22 3 

field. These provide a mean moment of 330 (+ 18)7 R ■ 4.8x10 gauss-cm 

at a tilt angle of 14°+5 and a longitude of 14S°+21°. This means 

tnac the surface field at Mercury is about 1*/. of Earth's, while tnc 
-4 

moment is 6x10 of Earth's. The polarity sense is the same. 

The origin of the field cannot be uniquely determined. It may 
be due to an active dynamo, a remanent magnetic field or a combination. 
Consideration of remanence as the source leads to some difficulties, 
although definitive knowledge of the planetary interior structure 
and thermal state is lacking sufficient to absolutely eliminate 
this source. Success in attempting to explain the field as due 

to an active dynamo has encouraged these efforts . Therefore, Mercury raa^ 
join Earth and Jupiter as an example of a planet possessing an 
rt. ..od fluid region with a convecting motion which regenerative.y 
maintains the magnetic field. The source of convective energy may be 
radiogenic decay and heat release, gravitational settling and 
differentiation or precessional torques. 
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INTRODUCTIO N 

Direct observations of the magnetic field of Mercury and its 
magnetosphere, as formed by the interaction of the solar wind, were 
made by the USA Mariner 10 spacecraft, in March of 1974 and again 
in March 1975 (Ness et al., 1974; Ness et al . , 1975; Ness et al., 
197o). From these data, it was discovered that there exists an 
intrinsic magnetic field of the planet Mercury, sufficiently strong 
at present to deflect the solar wind flow around the planet, thereby 
forming a well developed detached bow shock wave in the super- 
Alfvenic solar wind flow. Analyses of two other closely related 
experiments, the plasma and energetic particle instruments, have 
also appeared (Ogilvie et a., 1974; Hartle et al., 1975; Ogilvie 
et al., 1977 and Simpson et al., 1974), One unique aspect of the 
observations during the first encounter was the detection of a 
sudden change of the state of the magnetosphere, which was interpreted 
to he similar to that identified in the terrestrial magnetosphere 
as a 3ubstorm (Siacoe, Ness and Yeates, 1974). 

It is the purpose of this paper to very briefly review the 
magnetic field observations and their anal>',es as they relate to 
the determination of the magnetic field of the planet. In addition, 
the implication of this magnetic field with respect to its origin 
within the planetary interior is considered from two model viewpoints; 
passive remanent magnetization or an active core dynamo. The 
existence of a magnetic field at Mercury provides the only firm 
experimental evidence that Mercury is a differentiated planet, 
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consistent with theoretical models of planetary evolution which 
predict that to be the case. 

MAGNETIC FIKLU OBSERVATIONS 

Earlier publications have presented graphically the observations 
during the two encounters with Mercury by Mariner 10 as well as 
discussed their specific characteristics. Only a summary of salient 
features and some recent work extending the analyses is presented 
here. Figure 1 presents two views of the sub-Marlner 10 trace on 
the surface of the planet, while the spacecraft was located within 
its magnetosphere. Both close fly-bys occurred on the night-side 
of the planet and the third was especially selected to compliment 
the measurements and interpretations of data obtained from the 
first. The second encounter in September 1974 is not discussed in 
this paper because no direct observations of the magnetosphere or 
bow shock were possible since the spacecraft passed more than 
50,000 km from the surface of the planet, in order to provide good 
imaging coverage of the south polar region. 

Readily evident in Figure 1 is the contrast in the latitude 
region covered by the two fly-bys. The first was almost a pure 
equatorial pass, with a closest approach distance R * 1.30 
(planetocentric) while the second was almost a polar pass, with a 
closest approach distance of only 1.13 R», 

The data obtained from these two encounters is shown in 
Figure 2, where the magnitude of the field, F, and the relative 
fluctuations, 5/F, are presented on logarithmic scales. Individual 
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data points were obtained at i_'5 Hz and component averaged over 6 
second intervals, from which reconstructed vector field representations 
were obtained. Readily evident in this figure Is the significant 
difference in the magnitudes of the field measured in the first and 
third encounters, as well as the magnitude of the relative f luctuat ior s . 
The maximum field observed during the third encounter was 400 nT while 
that during the first encounter was only 100 nT, Itself a factor of 4-5 
times the interplanetary field strength. More important, notice 
should be made that following data point 90 of the first encounter, 
there is a significant increase in the relative fluctuations 
accompanied by a decrease in the magnetic field intensity. This 
signals the onset of the disturbance referred to earlier and already 
discussed as a Hermc-an substorm. 

An imporcant aspect to consider In the analysis and lnterpret.it ion 
of these data, with respect to the Intrinsic planetary magnetic field, 
is the rel itive contribution due to currents external to the planet 
which arise because of the deflected solar wind flow. Although 
specific measurements of the ion component of the solar wind plasma 
were not made by Mariner 10, measurements of the electron component 
of the plasma indicate differences in the momentum flux observed 
between the two encounters and the magnetic field observations also 
show differences in the interplanetary field. 

Thus, with different "upstream" solar wind conditions, it would 
not be expected that the resulting external current patterns would 
be identical for the two sets of data. This creates a problem in 
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tli.it one must choose which data set Is most representative of a 

Btat >nary magnetosphere . Clearly, the third encounter satisfies 

tliat requirement and carries with it the additional benefit that 

the average magnitudes measured during the third encounter were 

much larger than those during the first encounter. This is due partly 

to the fact that the average distance of third encounter measurements 

was 1.41 i< M as opposed to 1 76 for the first encounter. 

The distribution of field magnitude observed during the first 
and third encounters and the distribution of the radial distanacs 
of the spacecraft for the data sets I and III are shown in Figure 3. 
There is no question but that the two data sets are distincly 
diff< nnt in the magnitude distribution, where it Is seen that 
the median field magnitude lor the third encounter (1 d5 nT) is a 
factor of 2.5 times that of the median for the first encounter. 
INTRINSIC T1ELU ANALYSES 

Four different methods have been used to analyze the observations 
and derive quantitative values for the description of the planetary 
magnetic field; 

1. Comparison of bow .hock and magnetopause positions relative 
to those observed at Earthy 

2. Direct spherical harmonic analysis of the data, 

3. Modelling of the magnetosphere by use of an image dipole 
and infinite current sheet in addition to the planetary 
field, and 

4. Scaling of a mathematical model for the terrestrial 


magnetosphere 


f 
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Fhe results obtained yield dipole moments, M, from 2.4 to 

22 3 

5.1x10 gauss-cm , with equivalent equatorial Held Intensities of 
165 - 350 nT. In terms ot the individual gausslan coefficients 
the dipole moment Is defined as 

m - v^V + <*jV + a^V 

The lower values of the dipole moment are associated with those 
studies using methods 3 and 4 and employing partial quadrupole 
(g and octupole (g^ U ) terms to improve the least square fits 
of the model to the observations. In addition, users of methods 
3 and 4 have includ'd data from both encounter I ar.d III, thus 
mixing the data sets and the external current systems. 

Ibe bow shock am' magnetopause are located much closer to 
Mercury than at Earth, by a factor of 7.5, after normalizing by 
the planetary radius. This yields, along with the plasma observations 
of the momentum flux, an estimate of the equivalent dipole moment 
of 3-7x10 ^'gauss-cm by method 1. 

The data set avail able for use of method 2 is Incomplete, in 
the mathematical sense, because it docs not provide vector measure- 
ments over a closed surface sui rounding, the planet. Such a spatial 
distribution of data Is necessary in order to uniquely determine 
the separate contributions of magnetic fields due to sources 
internal to the planet and those external to the surface on which 
the measurements are performed. Thus, no unique determinations of 
the planetary magnetic field repn sentation can be derived rflth 


method 2 


Hu- use of only one of the five possible quadrupole moment 
terms and one of the eight possible octupolc moment terms in methods 
3 and 4 if. arbitrary and can not be justified. The use of these 
terms corresponds to a displacement of the dipole along its axis 
given by the following relationship: 



It is well known that n> ither the Karth nor Jupiter (Acuna and Naas, 
1976) show the axial synmetry of their planetary fields which would 
Justify incorporation of only the g ° and terms. 11>e ir inclusion, 
as used In methods 3 and 4 by Whang (1977), Jackson and beard (1977) 
and Ng and Beard (1978) yield offsets of the dipole, A^ ■ 0.2 K M , 
where Z is parallel o the dipole axis. 

A sunmary of all analyses of the planetary field at Mercury 
is presented in Table I. Included in the last entry are the results 
obtained in this paper and to be discussed shortly. Some of the 
differences in the dipole terms probably reflect the limitations 
of the models used to represent the external current systems and 
the confusion caused in the analyses by incorporating both data 
sets from the two encounters In addition, it is also certain that 
because of the relatively short radial excursion of the data, the 
putative g,° and gj 11 terms are also jpatially aliased w<th the g^ 0 
term. 

Recent analyses by method 2 of subsets of data from the third 
encounter taken near closest approach have shown a convergent series 
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of dipole moment values, which arc believed to best represent the 
Intrinsic planetary field. ihese results an shown In Figure 4, 
where a series of six successive solutions Is shown for each model, 
classified by 1m En. The value of m Indicates the highest order 
of the Internal terms whil_* the v iue of n Indicates the highest 
order of the external t rms. Note that as the subset identification 
number increases from 1 to 6, the number of data points Incorporated 
Increases from 20 to 125, It is seen that for all three solutions 
I1E0, 11 El and IlE^, using three different assumptions for the 
complexity of the external current system, the orientation of the dipole 
mcment vector Is well determined. The latitude Is found to be 
14°+5' while the longitude Is l4F n +21°. Hie error bars graphically 
Illustrate the orientation uncertainty In the main dipole moment 
term. Notice should bo taken of the stability of the equatorial field 
Intensity, which Is 332+4 nT lor all three solutions, while the root 
mean squared deviation of the fits of the models to the data 
decreases substantially from 19 to 9.3 nT as n Is increased. 

An Important consideration in th« analysis of such data by 
least squares approximation is whethci or not the analysis is well 
conditioned with respect to mathematical stability. A measure of 
this can be obtained by determining the condition number of the 
solution matrix. Figure 4 shows that the condition number 
Increases from 2 to >19 as the order of the external fipld Is increased. 


From experience, it is known that values of the condition number 
greater than 300 can lead to serious mathematical instabilities and 
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•o no higher order approximation* have been attempted. 

Hie data Intervals used In the Individual data sets are Illustrated 
In Figure 5, where the three components of the observations are 
shown compared to the IlKl model, previously published. Figure 5 
also includes a measure of the relative contribution of external source 
fields to the observations In the upper most panel on the left hand 
side. Ill ere it is seen that, for most of the data, the principle 
contribution to the observations is due to the Internal magnetic 
field. It is necessary to be cautious In any modelling process 
whereby one can improve a fit to a data set while invoking a physically 
Implausible model. Such a case could be developed here if one 
increased the order of the external fields beyond n ■ 2. 

As shown In Table I, those models which incorporate only a 
dipole term but use methods 3 and 4 provide dipole moments which 
are approximately 80/ of the value obtained by direct spherical 
harmonic analysis. Considering the Incorporation of the tuo encounter 
data sets and the differences caused by different solar wind conditions, 
thl9 is felt to be reasonably good agreement. However, when those 
methods incorporate higher order specially selected quadrupole and 
octupole terms, there is a substantial decrease in the dipole 
moment subsequently determined. As previously discussed, this is 
felt to be an artifact of the analyses and not to reflect the true 
character of the planetary field. 

INTER 10 d OF MERCURY 


A recent review of previous studies of the interior of Mercury 
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has bet’ll presented by Ness (1978). Here the results of those studies 
on the Internal structure * < summarized In Table II. The radius 
of the core formed during planetary evolution Is seen to be quite 
large, approximately 0.75 of the planet Itself, due to the large 
Iron nuiss fract Ion, deduced by all Investigators, t>57.. The possible 
internal structure of the planet, as summarized In Table II, Is 
based upon a number of assumptions tbout the conditions at the Lime 
of formation of the planet and Its sabsoquent evolution, The existence 
of a planetary magnetic field argues strongly for a differentiated 
planet since at present It is bellev.d that only an active dvnamo 
or a passive remanent magnetic field can he the source of any 
planetary field. Another common feature of the thermal evolution 
models is that the sub Curie poLnt depth is very close to thp 
planetary surface so that only a thin lithospheric shell is available, 
if remanence is to be the source of the external field. 

Table III suirsnarizes the studies conducted on remanent 

magnetization as a possible source of the planetary field of Mercury, 

3 22 3 

assuming the dipole moment of 330 nT R M »!».8xl0 gauss-cm . All 
the thermal evolution models have led tc lithospheric shells below 
the Curie point of only a few hundred kilometers thickness. Thus, 
the magnetization levels required, as shown in Table III, exceed 
sign! f leant! y those observed in the average lunar sample magnetization. 
Thus, remanent magnetization does not appear to be a very plausible 
source of the magnetic field ef Mercury. 

Therefore, the possibility that an a:tive dynamo must exist 
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at Mercury appear* compelling. To examine this, we now consider 
current developnuMit a In dynamo theory, also reviewed at this conference. 
The most promising model Which lias been developed Is that of Busse 
(1978), who Invokes contiguous, cylindrical convection cells parallel 
to the axis of rotation of the planet. He has demonstrated their 
existence In his analogue model studies and has also mathematically 
attempted ''heir analysis. Table IV sutmarizes the salient physical 
parameters for the three planets in the solar system known 
unequivocally to possess a magnetic field of Internal origin. The 
important point in this brief summary of Busse' s extensive work 
is that tin upper hound on the relative magnetic energy density factor, 
\/K*M, is approximately the same for all the three planets. Busse 
suggests, thusly, that his model has general applicability to 
the planets and the process ot dynamo formation. It is Inappropriate 
in this paper to review dynamo theory further and the reader is 
referenced to other articles in this conference proceedings for 
mor i irticlpmt* in that field. 

SUMMAKY 

As a result of unique observations by the USA Mariner 10 space- 
craft in 1974 and 1975, Mercury has been discovered to possess a 
modest but significant intrinsic planetary field. Because the field is 
weak, the solar wind greatly distorts the planetary magnetic field 
as it forms a magnetosphere. Quantitative analyses indicate that 
the undistorted plaretary field is primarily a dipole with it3 axis 
14°+5° from the rotation ax*s of the planet and in the same sense 
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as Earth's. The undistorted equatorial fle.d Intensity Is 330+18 nT 

22 3 

yielding an equivalent dipole magnetic moment of 4.8x10 gauss-cm . 

Tin* planet occupies an extremely large fraction of Its magneto- 
sphere, when compared to Earth. The average sub-solar point 
stagnation distance at the magnetopause is predicted to be 1.65+.15 Rj^, 
a factor of 6 times smaller than at Earth' s average position (10.4 Kg) . 

While thermo- remanent magnetization Is a possible source of 
the intrinsic field, It requires an anomalously high percentage of 
free-lton In the anticipated thin sub-Curie point lithosphere and an 
extremely strong an ent dynamo. Tills alternative appears somewhat 
implausible at present, and a more likely source Is an active dynamo. 

This requires a partially fluid core. Tills requirement Is consistent 
with theoretical models of planetary evolution which have been put 
forward. Indeed, the existence of a planetary f*eld argues 
convincingly that the planetary interior Is differentiated, regardless 
of the source of the magnetic field. 

Undoubtedly, remanent magnetic fields must exist on Mercury and 
when they are measured they may well provide an Invaluable historical 
record of the formation of the planet. Whether or not they shall be 
as revealing p s the studies of magnetic striping in the ocean bottoms, 
which led to the development of plate tectonic theory, is of course 
unknown at present. A very close orbiting spacecraft and/or a 
surface rover will be required In the future because the remanent 
fields may be so small. 
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1 NTEKIOR OK MERCURY 

KNOWN ESTIMATED 

P - 5.44 GMS/CC p (r , t ) T(r,t„) H(r,t ) 

o o o 

K - 2439 + 1 KM K(r,T) Cp(r f T) 

2 

t - 3 t - 58.65 DAYS f (DYNAMIC ELLI PTICITY) 

ROT ORB 


AUTHOR 

FE MASS 
FRACTION 

MF 

CORE 

R c /R 

C/MR 2 

PLAGEMANN (1965) 

- 

0.77 

0.86 

- 

KOZLOVSKAYA (1969) 

0.588 

0 

0 

0.368 


0 . 598 

0.598 

0.71 

0.324 

MAJEVA (1969) 

0. 58 

0 

0 

- 

REYNOLDS & SUMMERS 

U.709 

0 

0 


(1969) 

0 . 67 5-0 , 687 

0.681 

0.79 

- 

SIEGFRIED & SOLOMON 

0.69 

0 

0 

0.394 

(1974) 

0.62 

0 665 

0.75 

0.325 

TOKSOZ & JOHNSTONE 
(1975) 

— 

— 

0.69 

— 

SOLOMON (1976) 

0.62 

- 

0.75 

- 

FRICKER ET AL. (1976) 

0 . 66 


0.78 

tm 


TABLE II 
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TABLE III 


PROPERTIES OF PLANETARY CORES - BUSSE (1978) 



TABLE IV 


FIGUKF. CAPTIONS 


1. Two views of the trajectory of the Mariner 10 spacecraft as it 

transltted the llermean magnetosphere In March 1974 and March 1975. 
The planetocentrlc distance at entry and exit from the magneto- 
sphere Is shown, as is the point of closest approach, Identified 
by a + sign, 

iJ. Magnetic observations obtained by the Mariner 10 spacecraft 

during the two transits through the magnetosphere described in 
Figure 1. 

3. Statistical histograms of the field intensity and spatial position 
of the Mariner 10 spacecraft during the two encounters with the 
planet . 

4. Summary of spherical harmonic analyses to determine the magnetic 
dipole moment of Mercury. 

5. Comparison of observation and theory for analyses of the intrinsic 
planetary field during the Mariner 10 third encounter with Mercury. 
This figure Illustrates '.he relative position of the six subsets 
of data selected for sequential analysis. 
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